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ABSTRACT . . . . 

Microfracture (fiber/matrix loIdhT* compute- 

SSfStesSS ir 

that under thermal loads alone, microfracture propagation is not as sensitive as 
mechanical loads. 


INTRODUCTION 

High temperature metal matrix composite. (HTMMC) are 
for advanced propulsion systems, demanding high oper.t.on.1 tem^rature, ( 4 TOt. 1 00. C). 

High moduli .^strength, tailorable 

resistance make these materials especia y a , durability. Traditionally, researchers 

HTMMC is critical to assess their structural integrity and y- , 1 0 i 

"IkTat h m icmfracture using stresses, strain, and stress intensity f acto,, at the loca l level 

for the crack initiation and propagation. Such procedures are 

intensive and difficult to observe by conventional experiments. An alternate appro 
^r*e effect of microfracture on the global response (displacement, work done etc.). 

3L «UblLh the hierarchy of respective fracture modes for a 
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by fiber or matrix fracture. This demonstrates that debonding is a weaker fracture mode and 
very likely, it will instantaneously follow the stronger fracture modes (fiber/matrix fracture) 
when the composite is subjected to longitudinal tensile loads. Typical strain energy release rate 
curves for a unidirectional SiC/Til5 0.35 fvr composite subjected to longitudinal load are shown 
in figure 1. It shows that debonding mode of fracture requires less energy to propagate and 
occurs following either the fiber or the matrix fracture. Microfracture was also simulated for 
other types of mechanical loading. Similar observations were made for other types of mechanical 
loading in both unidirectional and crossply composite laminates. 

The objective of the present paper is to outline a procedure to computationally simulate 
microfracture and identify /quantify the microfracture modes and propagation for unidirectional 
and crossply metal matrix composite laminates subjected to thermal loads. 

FINITE ELEMENT MODEL 

The finite element models used in the computational simulation procedure consist of a 
group of nine fibers in a three-by-three unit cell array ("nine cell model"). The unidirectional 
composite system consists of 35 percent fiber volume ratio (fvr) SiC/Til5 metal matrix compos- 
ite (silicon carbide fiber and titanium alloy matrix). There are 16 nodal segments ("bays") 
along the length of the fiber. Each unit cell as shown in figure 2, consists of 40 hexahedron 
(six-sided) and 8 pentahedron (five-sided) solid elements for a total of 6912 elements and 6953 
nodes in the model. The crossply composite has three plies with 0/90/0 lay-up, as shown in 
figure 3, and is a 30 percent fvr SiC/Til5 metal matrix composite. There are six bays along the 
length of the fiber in both 0° and 90° plies. Each unit cell consists of 40 hexahedron (six-sided) 
and 8 pentahedron (five-sided) solid elements for a total of 2592 elements and 2863 nodes in the 
model. The properties of the constituents at the reference (room) temperature are shown in 
table I. 

In this composite, the coefficient of thermal expansion (CTE) of the matrix is more than 
that of the fiber (a m > a f ). The matrix longitudinal stress due to the cooldown from processing 
to use temperature is tensile. The corresponding fiber longitudinal stress is compressive. In the 
case of the crossply composite, there is also a mismatch between the thermal expansion coeffi- 
cient of the 0° and 90° plies, resulting in higher matrix microstresses. In any case, the matrix 
microstresses are tensile, while the fiber microstresses are compressive. Because of this stress 
state, the fracture is likely to initiate in the matrix. It is then propagated either through the 
fiber-matrix interface to cause debonding or through the interply layer to cause delamination. In 
the simulation, fracture is introduced around the fiber, such that the whole fiber circumference is 
debonded. Fracture is simulated by placing duplicate node points on either side of the crack. 
These duplicate nodal or grid points have the same geometrical location, but no connectivity 
exists between them, thus, in effect producing a crack of zero width. Symmetry boundary condi- 
tions are applied in the middle planes, so that the composite is free to move on either side when 
thermal loads are applied. Resulting stresses and internal strain energy corresponding to those 
applied thermal loads is computed by the finite element analysis. The corresponding strain 
energy release rates are computed for perturbed fracture configurations by using the definitions 
described below. Strain energy release rates (SERR) for different fracture modes are then 
compared to establish the fracture process, identify /quantify the various fracture modes and 
establish the hierarchy of those fracture modes. 
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STRAIN ENERGY RELEASE RATE 


The strain energy release rate (SERR) is an acceptable indicator of the fracture toughness 
of a material. It gives a measure of the amount of energy required to propagate a defect in a 
material. Hence, one can make a direct comparison of damage tolerances between different 
microfracture configurations (modes/paths), materials and geometries. One of the methods used 
to calculate SERR is the crack closure method. In this method, nodal displacements and 
corresponding nodal forces at the crack tip location, are used to determine the amount of work 
required to close the crack, which has been extended by an incremental amount. 

The above approach is a local or microfracture approach since the amount of energy 
produced by the local displacements and forces at the crack tip location are used ^to _calculate the 
corresponding strain energy release rates. Another approach to calculate global SERR is to 
compare internal strain energies of small fracture propagations. The strain energy release ra 

then, calculated as: 


dW _ 1 (S.E.); - (S.E.)! (i) 

dA 2 AA 


where 


dW incremental work done 

A A area of the new surfaces generated 

(S.E.)j, (S.E.) 2 strain energy prior to and after A A, respectively 

The SERR were computed by using both the crack closure method and by using global 
strain energy for one case to compare the values. Both methods give the same results, although 
using the global strain energy approach is computationally more effective and direct. However, 
the crack closure method is required to identify each failure mode in the presence of combined 
mode fracture. The global strain energy approach is used in the present work. The advantage 
of using the global SERR approach is that it bypasses local stress details to describe stress 
gradients that usually require relatively fine meshes. 


MULTIFACTOR INTERACTION RELATIONSHIP 


Constituent material properties are required at higher use temperatures. Constituen 
material properties at reference (room) temperatures are known, as shown m table I. Ihe ettect 
of temperature on material properties is taken into account by using a multifactor interaction 
relationship (refs. 3 and 4) as shown in figure 4. It models the material behavior using a 
time-temperature-stress dependence of the constituent’s properties in a "material behavior 

space," as follows: 



( 2 ) 
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where 


P property 

T temperature 
S strength 

o stress 

0 reference 

F final 

m,n are exponents 

It assumes that various factors such as temperature, stress, stress rate, etc. influence the in-situ 
constituent material behavior. The multifactor interaction relationship (eq. (2)) represents 
gradual effects during most of the range and rapidly degrading properties near the final stages as 
has been observed experimentally (ref. 5). The exponents are determined from experimental 
data, wherever possible, otherwise default values for exponents are used which are based on 
"best judgment" from studies conducted on other materials (ref. 5). 

In the present investigation, in-situ properties are assumed to depend only on temperature 
(m = 0). The value of exponent n is taken as 0.5 for matrix and 0.25 for fiber. The final 
temperature is assumed to be the melting temperature of the constituent and the reference 
temperature is taken as the room temperature (21 °C/70 °F). For example, to calculate a 
matrix property at 815 °C, where the titanium alloy matrix melting point is taken as 927 °C 
(1700 °F), then from equation (2): 


P 

Po 


10.S 


927 - 815 
927 - 21 


0.35 


Any matrix property, i.e., modulus, at 815 °C will be 35 percent of the reference property, 
i.e., room temperature modulus. In the case of a, coefficient of thermal expansion, the exponent 
n is assumed to be negative, because o is assumed to increase with increase in temperature. 
Constituent material properties at 815 °C (1500 °F), calculated using equation (2) are shown in 
table II. 


RESULTS/DISCUSSION 

In the first set of simulations, both unidirectional and crossply composites were uniformly 
heated from room temperature (21 °C/70 °F) to a temperature of 300 °C (570 °F), i.e., AT of 
280 °C (500 °F). Constituent properties at room temperature were used and were assumed to 
remain constant for this thermal loading case. In the case of the unidirectional composite, 
fracture was initiated in the matrix at the middle plane of the center cell and then propagated 
through the fiber-matrix interface. The change in total strain energy was very small. Even 
when all the fibers were debonded, the change in total strain energy and the corresponding 
SERR were negligibly small (G = 0.02 lb/in.). In the case of crossply composite, fracture was 
again initiated in the matrix at the middle plane of the center cell and propagated along the 
fiber-matrix interface. When the center cell fiber was fully debonded, the change in total strain 
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energy from reference (no fracture) state to this fracture mode, was 0.4 percent. The corre- 
sponding strain energy release rate was also negligibly small. Fracture initiated at that location 
wm also propagated to the interply layer in order to simulate delamination between the top an 
middle pfies When the delamination extended up to four bays through the width (y-direction) 
and in the x-direction, there was no noticeable change in the total strain energy as comp " 
the reference state. Hence, it can be concluded for both unidirectional and crossply c °mp<*“tes, 
that thermally driven microfracture propagation is quite insensitive to temperature increases up 
to 260 ®C (500 °F) from the reference (room) temperature. 

In the next set of simulntions, the composites were cooled down from 815 “C (1500 ®F) to 
185 “C f 300 "FI ie a AT of -1000 “C (-1800 °F). Constituent material properties were 
computed at 815 "C (1500 “F), shown in table 11, by using equation (2) and were aammed to 
remL constant through this thermal hmd. Since, the CTE increas* ■ -■* * 
tlire the mismatch in CTE of the fiber and the matrix is very high at 815 C (table 11). As 
mentioned before, the thermal expansion coefficient of the titanium alloy (Til5) matrix is more 
than the thermal expansion coefficient of the silicon carbide ( S|C ) Fiber. us, UM “® * 
down process from high temperatures to cryogenic temperatures, the fiber eve ops p 
longitudinal stress, while the matrix is under tensile longitudinal stress. Therefore, the 
microfracture is likely to initiate in the matrix. In the case of the unidirectional comfmsite the 
microfracture was initiated in the middle plane of the center cell matrix and propagated n the 
matrix or through the fiber-matrix interface. The SERR curve for 

the fiber-matrix interface is shown in figure 5. Hence, it can be concluded that once approxi- 
mately 10 percent of the fiber surface is debonded, debonding can propagate at the same energy 
level. The SERR for the case when the fracture propagates in the matrix is zero. 1 hus, the 
fiber-matrix interface debonding is the only likely mode of fracture propagation in this case 
The unidirectional composite, therefore, will exhibit matrix cracking or a significant amount of 
interface debonding and is unlikely to fail in a brittle manner. 

In the case of the crossply composite, cracks can initiate in the 90° ply as intraply cracks 
or through the thickness of the 90° ply as transply cracks. When the crack initiates within he 
90° as an intraply crack at the mid-length of the fiber (fig. 6), it propagates in the matrix and 
then the fiber-matrix interface to debond the adjacent center cell fiber. T e curve as 

center cell fiber debonds is shown in figure 6. The curve shows the same pattern as observed 
earlier, once about 10 percent of the fiber surface is debonded, debonding can P ro P a 8 ate a 
same energy level. If the crack initiates in the matrix as a transply or through-the-thickness 
crack at the center of the specimen, then it can propagate through the matrix to debond the 
fiber in the adjacent (0°) ply, or it can propagate through the interply layer to cause termina- 
tion. When the transply crack propagates in the interply layer, there is no change m g 
strain energy and thus, the SERR is also zero. However, if the transply crack propagates 
through the matrix and debonds the adjacent ply fiber, the SERR for that fracture mode is 
shown in figure 7. As before, the fiber-matrix interface debonding is the only mode of fracture 
propagation in this composite. Hence, it can be concluded that thermally-driven microfracture 
in crossply composites will exhibit matrix cracking and fiber debonding due to the cool down 
from high temperatures to cryogenic temperature and the microfracture will not propagate in the 
interply layer to cause delamination. 

Other thermal loading cases involving thermal gradients were also evaluated, including one 
in which the top surface was maintained at a very high (815 “C/1500 "F) temperature and 
bottom surface was maintained at cryogenic (-300 "F/-185 “C) temperature. In other loading 
case top and bottom surfaces were maintained at 538 “C/1000 “F mid the middle plane was 
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kept at 815 °C/1500 °F. Both of these cases did not exhibit any sensitivity for microfracture 
propagation in either unidirectional or crossply composite. The results for microfracture ini- 
tiation and propagation are summarized in table III. For example, in the case of the unidirec- 
tional composite subjected to a uniform temperature increase, the microfracture is likely to 
initiate in the matrix and propagates through the fiber-matrix interface. If the microfracture 
initiates in the interphase region or the interply region (pre-existing crack), it does not propa- 
gate when subjected to such a thermal load. 


CONCLUSIONS 

A computational simulation procedure is described for microfracture initiation and 
propagation in metal matrix composites subjected to thermal loading. The procedure is applied 
to evaluate microfracture for both unidirectional and crossply composites. The significant 
results from this investigation are as follows: 

1. The crack closure method and the global strain energy methods give the same results for 
strain energy release rates (SERR). The global strain energy method is computationally more 
effective. However, crack closure method should be used to identify the individual modes of 
fracture in cases of mixed mode fracture. 

2. In general, microfracture propagation in HTMMC is not as sensitive to thermal loads as 
it is to mechanical loads. 

3. For both unidirectional and crossply composites, microfracture propagation is not 
sensitive for a temperature increase of 260 0 C/500 °F from room temperature. 

4. When cooled down from high temperature to cryogenic temperature, the fracture 
initiates in the matrix: 

i. For unidirectional composites, fiber-matrix interface is the only likely mode of fracture 
propagation. 

ii. For crossply composites: 

a. If the crack initiates within a ply between the fibers (intraply cracks), it propagates 
to debond the adjacent fiber. 

b. If the crack initiates through the thickness of the ply (transply crack), it propagates 
to the adjacent ply and debouds the fiber in that ply. 

c. If the transply crack initiates closer to a fiber, it debonds the adjacent fiber within 
the ply. 
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TABLE I. - PROPERTIES OF CONSTITUENT 
MATERIALS OF SiC/Til5 AT 
ROOM TEMPERATURE 



SiC 

fiber 

Til5 

matrix 

Inter- 

phase 

Modulus, E (Mpsi) 

62.0 

12.3 

12.3 

Poisson’s ratio, v 

.3 

.32 

.32 

Shear modulus, G (Mpsi) 

23.8 

4.8 

4.8 

Coefficient of thermal 
expansion, a, (ppm/°F) 

1.8 

4.5 

4.5 


TABLE II. - PROPERTIES OF CONSTITUENT 
MATERIALS OF SiC/Til5 AT 
815 »C (1500 °F) 



SiC 

fiber 

Til5 

matrix 

Inter- 

phase 

Modulus, E (Mpsi) 

57.0 

4.3 

4.3 

Poisson’s ratio, v 

.28 

.15 

.15 

Shear modulus, G (Mpsi) 

22.4 

1.9 

1.9 

Coefficient of thermal 
expansion, a, (ppm/°F) 

1.96 

12.8 

12.8 
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TABLE III. - MICROFRACTURE INITIATION AND MODE HIERARCHY 

UNDER THERMAL LOADS 


Constituent 

Composite/Load ' 

Fi 

ber 

:o 

Matrix 

(m) 

Interphase 

0) 

Interpl 

0 

y layer 
p) 

l a 

El 

I 


I 

n 

I 

p 

Unidirectional 

composite 

(UDC) 

Uniform 

temperature 

■ 

■ 

m 

■ 


m 

»P 


Thermal gradient 
through-the- 
thickn ess 

1 

1 

m 

l 



ip 


Crossply 

composite 

(CPC) 

Uniform 

temperature 

■ 

■ 

m 

■ 

i 

m 

ip 


Thermal gradient 

through-the- 

thickness 

1 

1 

m 


i 


ip 



*1: microfracture initiation. 
b P: microfracture propagation. 
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Figure 3. — Schematic diagram of crossply composite nine cell 
model. 
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Rationale: 

• Gradual effects during most range, rapidly degrading near final stages 

• Representative of the in situ behavior for fiber, matrix, interphase, coating 

• introduction of primitive variables (PV) 

• Consistent in situ representation of all constituent properties In terms of PV 

• Room-temperature values for reference properties 

• Continuous interphase growth 

• Simultaneous interaction of all primitive variables 

• Adaptability to new materials 

• Amenable to verification Inclusive of all properties 

• Readily adaptable to incremental computational simulation 


Subregions of 
intralaminar 
nonunfformfty — ( 



Notations: 

P - property; T - temperature; S - strength; R - metallurgical reaction; N - number of cycles: 
t - time; over dot - rate; subscripts: 0 - reference; F - final; M - mechanical; T - thermal 

Figure 4.— Assumed multifactor Interaction relationship. 
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Figure 5. — Strain energy release rate for the unidirectional com- 
posite; cooling from 1500 °F to -300 °F; matrix Initiated crack. 
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Figure 6. — Strain energy release rate versus fiber surface de- 
bonded for crossply composite, cooling from 1500 °F to 
-300 °F; intraply crack. 


6th Bay —v 



20 40 60 60 100 
Fiber surface d ©bonded, percent 


Figure 7. — Strain energy release rate versus fiber surface debond- 
ed; crossply composite under uniform cooling from 1500 °F to 
-300 C F; transply crack. 
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